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It  is  difficult  to  model  the  behavior  of  the  battery  pack  accurately  due  to  the  electrochemical  charac¬ 
teristics  variations  among  cells  of  a  battery  pack.  As  a  result,  accurate  state-of-charge  (SOC)  and  state-of- 
health  (SOH)  estimation  for  the  battery  pack  is  a  case  provocation.  The  estimation  process  poses  more 
challenges  after  substantial  battery  aging.  This  paper  tries  to  estimate  the  SOC  of  a  Li-ion  battery  pack  for 
an  electrical  vehicle  using  improved  extended  Kalman  filter  (IEKF)  which  benefits  from  considering  aging 
phenomenon  in  the  electrical  model  of  cells.  In  order  to  assemble  a  battery  pack,  we  find  cells  with 
similar  electrochemical  characteristics.  Model  adaptive  algorithm  is  applied  on  the  corresponding  cells  of 
a  string  to  minimize  cell-to-cell  variation’s  effect.  During  the  operation,  the  values  of  electrical  model  of 
each  cell  are  updated  by  the  same  algorithm  to  compensate  aging  effects  on  SOC  estimation  error.  The 
mean  value  of  updated  cell’s  model  is  used  for  a  single  unit  cell  model  of  the  pack  used  at  IEKF  to  achieve 
more  accurate  SOC  estimation.  The  algorithm’s  fast  response  and  low  computational  burden,  makes  on¬ 
board  estimation  practical.  The  experimental  results  reveal  that  the  proposed  approach’s  SOC  and  voltage 
estimation  errors  do  not  exceed  1.5%. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Large-scale  energy  storage  consisting  rechargeable  batteries, 
have  become  one  of  the  most  popular  candidates  for  electrical 
energy  storage  for  several  reasons.  Matching  supply  and  demand  in 
power  grids,  providing  fast  response  to  energy  demand,  ease  for 
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siting  and  their  high  energy  efficiency  are  some  to  name.  Under  the 
global  demand  for  reduction  in  greenhouse  gas  emissions  in  the 
power  sector  as  well  as  car  industry  and  emergence  of  renewable 
energy  sources,  advanced  battery  systems  are  proposed  for  wide 
range  of  applications  varying  from  electrical  vehicles  (EVs)  and 
hybrid  electrical  vehicles  (HEVs)  all  the  way  up  to  smart  grids. 
Nowadays,  Li-ion  batteries  have  become  the  preferred  choice  in  the 
field  of  power  battery  research  and  EVs  applications.  They  have 
gained  this  popularity  as  a  research  topic  in  industry  and  academia 
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because  of  the  characteristics  of  high  energy  density,  less  weight, 
longer  cycle  life  and  no  memory  effect  [1  .  Storage  systems  for  EVs/ 
HEVs  applications  to  attain  their  voltage  and  traction  power  re¬ 
quirements  include  tens  to  hundreds  of  Li-ion  single  cells  con¬ 
nected  in  series  and  parallel  or  series/parallel.  Moreover,  to  have  a 
secured  and  safe  energy  source  and  distribute  power  in  EVs/HEVs  a 
battery  management  system  (BMS)  needs  to  have  an  accurate  inline 
estimation  of  SOC  and  SOH  for  Li-ion  cells  in  the  battery  pack  [2]. 
An  appropriate  battery  model  is  necessary  for  proper  design  and 
operation  of  battery  systems  using  BMS.  This  battery  model  also 
helps  facilitating  estimation  schemes  stated  above  [3  .  Several 
modeling  approaches  are  available:  empirical  models,  statistical 
models,  and  circuit  models  [4-8  .  Since  applications  covered  in  this 
paper  are  related  to  EVs/HEVs  and  smart  grid,  an  electrical  model  is 
preferable  [9  .  There  are  numerous  models  proposed  in  the  litera¬ 
ture  which  are  accurate  enough  for  showing  the  electrical  behavior 
of  Li-ion  batteries  [9-12  .  Yet,  these  models  rely  on  internal  pa¬ 
rameters  of  the  battery  such  as  the  SOC,  which  is  hard  to  measure 
[13].  Accurate  estimation  of  SOC  under  any  operating  condition  still 
remains  a  challenge.  The  estimation  process  even  poses  more  dif¬ 
ficulty  in  the  battery  pack  due  to  cell-to-cell  variations. 

In  recent  years,  there  has  been  great  effort  for  improving  the 
accuracy  of  SOC  estimation  of  single  cells  and  battery  packs  [14  .In 
the  case  of  single  cells,  the  coulomb  counting  method  is  the  most 
common  method  used  for  SOC  estimation  [15  .  However,  this 
method  has  several  drawbacks  including  sensitivity  to  the  initial 
SOC  value  which  leads  to  inaccurate  estimation  and  accumulated 
error  caused  by  use  of  integration.  There  are  several  new  and 
improved  methods  which  are  based  on  artificial  neural  networks 
(ANNs),  fuzzy  logic,  adaptive  observers,  sliding  mode  observers, 
extended  Kalman  filter  (EKF)  and  unscented  Kalman  filter  (UKF) 
[13,16-20]  or  integration  of  several  individual  approaches  like 
coulomb  counting  and  EKF  [21  .  Last  estimation  methods  are  im¬ 
provements  to  classical  Kalman  filter  (KF)  which  are  suitable  for 
nonlinear  systems  18-21  .  Among  the  SOC  estimation  methods, 
some  are  used  in  BMSs  in  EVs/HEVs  successfully  which  are  based  on 
EKF  [18].  However,  there  are  issues  regarding  accuracy  of  SOC 
estimation  as  a  battery  pack  ages  which  is  discussed  here. 

Batteries  lose  a  portion  of  their  capacity  in  the  process  of  aging 
which  also  causes  changes  in  the  battery  model  parameters  [22,23]. 
Moreover,  it  is  important  to  recognize  that  SOC  is  a  short  term 
characteristic  while  capacity  fade  is  a  long  term  one.  Yet,  they  have 
a  close  correlation  with  one  another.  In  order  to  correlate  the  SOC 
with  capacity  fade  and  changed  model  parameters,  a  method  for 
accurate  estimation  of  SOC  is  required.  Therefore,  SOC  estimation 
method  for  precise  prediction  needs  to  update  its  model  parame¬ 
ters  adaptively  during  aging.  Most  recent,  a  novel  algorithm  was 
proposed  to  update  the  parameters  of  a  Li-ion  cell  model  using  EKF 


in  which  highest  priority  is  given  to  estimate  each  cell’s  voltage 
accurately  [24].  An  adaptive  extended  Kalman  filter  (AEKF)  is  pro¬ 
posed  to  estimate  the  SOC  based  on  cell  voltage  measurement  [25]. 
However,  the  algorithm  focuses  on  covariance  matching  and  ig¬ 
nores  aging  phenomenon  in  cells.  Estimations  which  are  based  on 
unscented  Kalman  filter  (UKF)  or  sigma  point  Kalman  filter  (SPKF) 
[19,20]  are  estimating  physical  model’s  parameters  as  well  as  SOC, 
however  this  convergence  takes  a  considerable  amount  of  time  to 
reach  its  final  value. 

As  the  above  methods  do  not  consider  cell  differences,  they  are 
not  applicable  for  the  battery  pack  and  mostly  focus  on  single  cell’s 
SOC  and  terminal  voltage  estimation. 

In  our  previous  work  [26  ,  brute-force  search  method  is  used  as 
an  optimization  algorithm  for  finding  new  electrical  model  for  an 
aged  cell  used  in  EKF.  Using  this  method,  SOC  estimation  error  is 
less  than  4%.  In  Refs.  27],  the  SOC  estimation  sensitivity  to  resistors 
and  capacitors  for  the  model  in  Fig.  1  is  investigated.  Based  on 
above  analysis,  it  is  clear  why  we  choose  resistors  in  voltage 
response  circuit  side  in  the  model  which  is  going  to  be  updated.  For 
the  vehicular  operation,  the  electrical  energy  storage  system  is 
composed  of  hundreds  of  cells  connected  in  series/parallel  due  to 
high  energy/power  demand.  In  this  paper,  a  string  cells  in  series  is 
considered  as  a  pack.  Concerning  the  battery  pack,  the  charging/ 
discharging  capability  of  weakest  cell  is  the  limiting  factor  in  the 
branch  [28].  Therefore,  the  characteristic  of  this  cell  for  battery  pack 
is  important  during  operation.  That  is  why  BMS  needs  to  be  aware 
of  each  cell’s  SOC  to  avoid  overcharging/over  discharging.  One  of 
the  stated  methods  is  screening  process  of  single  cells  and  use  of 
very  similar  cells  for  making  a  pack  [29].  However,  this  method  is 
time  consuming  and  lots  of  cells  would  be  omitted,  while  aging 
process  is  not  considered.  In  Ref.  [30]  an  estimator  for  Li-ion  bat¬ 
teries  is  used  which  estimates  SOC  for  each  single  cell  in  a  pack  of  N 
cells  connected  in  series.  The  method’s  computational  burden  for  a 
pack  of  cells  is  heavy  and  it  is  impractical  to  implement  it  for  in-line 
estimation  in  a  car  31]. 

This  paper  deals  with  design  of  battery  pack  which  is  composed 
of  120  cells,  connected  in  series,  providing  a  high  voltage  level.  A 
single  unit  cell  model  of  a  pack  is  employed  in  order  to  find  an 
accurate  SOC  estimation  for  both  new  and  aged  battery  pack.  The 
single  unit  cell  model  is  found  based  on  the  mean  values  of  updated 
cells’  model.  For  choosing  cells  with  similar  electrochemical  prop¬ 
erties  working  as  a  filter,  model  adapting  (MA)  method  approach 
introduced  in  Ref.  26]  is  used  in  which  filtered  cells  are  connected 
in  series.  MA  filtering  method  guarantees  cells’  similarity  in  the 
pack  and  consequently  enables  modeling  the  pack  with  a  single 
unit  cell.  Also,  the  same  approach  is  used  to  update  the  single  unit 
cell’s  model  as  the  battery  pack  ages.  Moreover,  the  improved 
extended  Kalman  filter  (IEKF)  which  improves  the  prediction 
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Fig.  1.  Schematic  diagram  of  the  electrical  model  based  on  [5]. 
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accuracy  by  changing  Kalman  gain  is  applied  to  estimate  the  bat¬ 
tery  pack’s  SOC.  The  obtained  results  by  FT75  and  NEDC  reveal  that 
the  proposed  approach’s  SOC  and  voltage  estimation  error  for  an 
aged  pack  do  not  exceed  1.5%. 

The  rest  of  this  paper  is  organized  as  follows.  Section  2  presents 
an  electrical  model  used  for  LiFeP04  cell,  which  is  also  used  as  single 
unit  model  for  the  pack  of  Li-ion  cells  in  the  pack.  A  review  of  the  EKF 
is  presented  in  Section  3.  Model  adaptive  approach  which  is  used  for 
filtering  cells  with  similar  electrochemical  characteristics  is  dis¬ 
cussed  in  Section  4.  In  Section  5,  cells’  model  connected  in  series  for 
battery  pack  is  built  based  on  electrical  model.  The  IEKF-based  SOC 
estimator  is  proposed  for  SOC  estimation  of  the  battery  pack  using 
the  single  unit  cell  model.  The  very  next  section  describes  the 
experimental  test  applied  on  battery  pack.  In  this  study,  a  test  on  a 
pack  of  Li-ion  cells  which  contains  120  LiFeP04  cells  connected  in 
series  is  carried  out  to  verify  the  estimation  accuracy  of  model 
adaptive-IEKF  method  for  SOC  and  terminal  voltage  estimation.  The 
paper  concludes  with  a  summery  and  final  remarks. 

2.  Battery  electrical  modeling 

SOC,  as  one  of  the  most  important  information  in  BMS,  is  an 
inner  state  of  each  cell  [31  which  cannot  be  measured  directly 
during  battery  operation.  As  a  result,  estimating  SOC  is  the  only  way 
to  obtain  its  value.  To  estimate  this  value,  a  model  for  the  battery  is 
needed.  A  variety  of  battery  models  are  developed  to  capture  Li-ion 
battery  performance  for  various  purposes,  among  them  the 
equivalent  circuit  models  and  the  electrochemical  models,  are 
widely  used  in  EV  studies.  The  electrical  circuit  models  use  equiv¬ 
alent  electrical  circuits  to  show  I—V  characteristics  of  batteries  by 
using  voltage  and  current  sources,  capacitors  and  resistors.  Due  to 
the  remarkable  relaxation  effect  of  the  Li-ion  battery,  the  model 
requirements,  including  enough  accuracy  and  covering  different 
empirical  conditions  like  working  condition  of  EVs/HEVs,  the 
model  presented  in  Ref.  [9  is  selected  as  the  battery  model-shown 
in  Fig.  1.  Energy  balance  circuit  is  a  part  of  model  which  delivers 
SOC  to  the  voltage  response  circuit.  In  this  model,  the  ohmic 
resistance  R0  consists  of  the  bulk  resistance  and  surface  layer 
impedance,  accounting  for  the  electric  conductivity  of  the  elec¬ 
trolyte,  separator  and  electrodes.  The  activation  polarization  is 
modeled  by  Rs  and  Cs,  and  the  concentration  polarization  is  pre¬ 
sented  by  R\  and  Q. 

To  cover  all  practical  conditions  and  considering  suitable 
complexity,  model’s  components  are  assumed  to  be  function  of 
SOC,  C-rate  and  temperature.  However,  since  optimization  algo¬ 
rithm  used  in  this  paper,  inherently  has  some  errors  to  estimate 
model  parameters  with  very  high  accuracy,  discharging  equations 
presented  in  Ref.  [9]  are  neglected.  It  is  assumed  that  there  is  just 
one  operating  function  for  charging  and  discharging. 

The  electrical  behavior  of  the  practical  model  in  Fig.  1  can  be 
expressed  as  follows: 

Vt  —  ^oc  —  ^trans  —  (1) 


and  long  time  transient  voltage  responses  for  charging/discharging 
respectively. 

3.  Extended  Kalman  filter 

Kalman  filter  (KF)  is  a  well-known  estimation  theory  introduced 
in  1960  [32].  The  filter  provides  a  recursive  solution  through  a 
linear  optimal  filtering  for  estimating  systems’  state  variables. 
However,  if  the  system  is  nonlinear,  a  linearization  process  is  un¬ 
dertaken  at  each  step  which  is  used  to  approximate  the  nonlinear 
system  with  a  linear  time  varying  (LTV)  system.  Using  the  LTV 
system  in  KF,  would  lead  to  an  extended  Kalman  filter  (EKF)  on  a 
real  nonlinear  system  18].  The  calculation  process  for  a  nonlinear 
system  including  its  modeling  is  given  in  equations  (5)  and  (6): 

Xfc+1  =f(xk,uk)+wk  (5) 

yk+i  =  g(Xk,  Ufc)  +  vk  (6) 

where  equation  (5)  represents  all  of  the  system  dynamics 
expressed  in  state  equations  and  equation  (6)  indicates  the  output 
equation  of  the  system  with  a  static  relationship.  Function  f(xk,Uk)  is 
a  nonlinear  transition  function  and  g(xk,Uk)  is  a  nonlinear  mea¬ 
surement  function.  Vectors  wk  and  vk  denote  process  and  mea¬ 
surement  noise  which  are  uncorrelated  zero-mean  white  Gaussian 
stochastic  processes  with  covariance  matrixes  Qand  R  respectively. 

At  each  time  step,  matrices  of  f[xklUk)  and  g{xk,Uk)  are  linearized 
close  to  the  operation  point  by  the  first  order  in  Taylor-series  and 
the  rest  of  series  are  truncated.  Assuming  that  f{xk,Uk)  and  g{xk,Uk) 
are  differentiable  at  all  operating  points  and  Ak  =  df/dx |  -, 
Hk  =  dg/dx |  -;  EKF  starts  filtering  with  the  best  available  infor¬ 
mation  on  the  initial  state  (xj )  and  error  covariance  (P+)  as  shown 
in  equation  (7). 

Xq  =  E[x0],Pj  =  £[(x-Xo)(x-Xo)r]  (7) 

An  illustration  of  the  EKF  is  presented  in  Fig.  2  [33,34]. 

4.  Model  adaptive  approach  for  SOC  estimation 

For  the  efficient  management  and  control  of  battery  pack,  the 
BMS  needs  to  have  accurate  estimation  of  SOC.  EKF  is  an  optimum 
state  estimator  for  nonlinear  system  which  works  by  recursion.  EKF 
can  be  used  for  SOC  estimation,  since  EKF  is  formed  in  discrete-time 
state  space,  equations  (5)  and  (6)  should  be  transformed  into 
discrete  form.  Following  the  form  in  EKF,  the  state  equations  for  the 
nonlinear  system  of  the  battery  are  x\  =  Vs,  X2  =  V/,  X3  =  SOC.  SOC  in 
contrast  to  battery’s  terminal  voltage  Vt  cannot  be  measured 
directly.  Discrete  time  state  space  form  for  practical  model  after 
linearization  of  equations  (5)  and  (6)  are  shown  in  equations  (8) 
and  (9). 

Xfc+1  =  +  BklU{  +  wk  (8) 


Vtrans  —  Vs  +  V[ 


(2)  Vt  —  yk  —  Hkxk  +  DkILk  +  vk 


(9) 
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The  state  vector  of  practical  model  consists  of  three  state  vari- 
(3)  ables  as  shown  in  equation  (10). 


(4) 


Xk 


Vs,k 

Vl,k 

sock 


(10) 


where  Vt  is  the  battery  terminal  voltage,  Voc  is  the  battery  open  where  SOC/,.  is  the  observation  of  SOC  at  time  step  k  which  is  equal 

circuit  voltage  (OCV)  and  k  is  load  current.  Vs  and  V/  are  the  short  to  equation  (11): 
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Initial  estimates  for  xk_1  and  Pk_1 


Fig.  2.  Complete  picture  of  the  operation  of  the  extended  Kalman  filter. 
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where  At  equals  to  sampling  time,  is  the  columbic  efficiency  and 
Cusabie  is  the  battery’s  available  usable  capacity.  The  observation 
equations  of  the  discrete  system  are  as  follows: 
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4.2.  Model  adaptive  algorithm 

The  model  adaptive-IEKF  is  a  novel  approach  presented  in  this 
paper  for  battery  pack  SOC  estimation  and  it  is  usable  for  single 
cells,  as  well.  This  approach  consists  of  two  parts,  model  adaptive 
part  and  IEKF.  As  battery  ages,  estimated  SOC  by  EKF  loses  its  ac¬ 
curacy  [27].  For  precise  SOC  estimation  of  an  aged  cell,  EKF  requires 
an  updated  electrical  battery  model.  Sensitivity  analysis  performed 
in  Ref.  [27  shows  that,  R0 ,  Rs  and  R\  in  voltage  response  circuit  side 
of  the  model  have  the  most  effect  on  SOC  estimation  by  the  EKF.  In 
other  words,  as  Li-ion  cell  ages,  those  parameters  must  be  updated 
in  the  battery  model.  Fig.  3  represents  SOC-Vt  graph  for  a  new 
LiFeP04  cell  and  also  the  same  cell  after  aging,  discharged  with 
constant  current  (1C  in  this  case).  Mean  values  for  electrical  ele¬ 
ments  of  the  new  and  degraded  cell  are  presented  in  Table  1.  It  is 
obvious  from  Fig.  3  that  terminal  voltage  has  sharp  drops  and  rises 
at  specific  states  of  SOC.  In  the  case  of  examined  electrochemistry 
cell,  the  sudden  voltage  variations  happen  at  SOC  values  between 
15%  and  92%  which  are  considered  as  reference  SOCs  here. 

Assigning  already  found  two  values  of  SOCs  to  voltage  deriva¬ 
tive,  the  model  adaptive  approach  brings  further  innovation  in  the 
algorithm.  The  idea  is  implemented  in  replacing  the  obtained  new 
model  in  EKF  with  the  old  model. 

According  to  equation  (1 ),  Vt  is  the  open  circuit  voltage  (OCV)  of 
cell  deducted  by  Vtrans  and  voltage  drop  on  cell’s  internal  resistance. 
Since  parameters  R0l  Rs ,  Cs,  R\  and  Q  depend  on  cell’s  SOC,  Vt  is  a 
function  of  SOC  as  shown  in  equation  (13). 


SOC 

Fig.  3.  SOC-Vt  graph  for  a  new  LiFeP04  cell  and  the  same  cell  after  aging  (discharged 
with  1C). 


Vt  =  h(  SOC)  (13) 

Taking  derivative  of  Vt  in  equation  (13)  gives: 

Vt  =  SOC*/i(SOC)  (14) 

Knowing  that  SOC  =  — /L/Cusabie,  equation  (14)  is  reduced  to: 

Vt  =  *h(S  OC)  (15) 

^usable 

Graph  of  -Vt  for  a  new  cell  and  an  aged  cell  are  shown  in  Fig.  4. 
There  are  two  abrupt  changes  in  Vt  versus  SOC  graph  discharged 
with  a  constant  current  (1C).  These  two  points  are  assumed  refer¬ 
ence  points  in  which  their  SOC  can  be  determined  from  cell’s 
chemistry.  The  SOC  values  for  the  reference  points  are  set  to  92% 
and  15%. 

According  to  [27],  the  battery  model  is  sensitive  to  R0l  Rs  and  R\ 
in  electrical  circuit.  In  addition,  each  of  the  elements  are  a  function 
of  SOC  according  to  the  following  equations  [9]: 

R0  =  f^SOC4  +  b2  SOC3  +  b3  SOC2  +  b4SOC  +  b5  (17) 

Rs  =  c1e_C2SOC  +  c3  (18) 

R,  =  g^soc  +  g3  +  g4SOC  (19) 

In  equations  (17)— (19),  65,  C3  and  g3  have  the  most  effect  on  R0 , 
Rs  and  Ri  values.  Assuming  all  parameters  to  be  constant  except  for 
fi5,  c3  and  g3  in  the  aging  process  of  battery,  the  variable  ones  are 
updated  using  proposed  optimization  algorithm. 


372 


S.  Sepasi  et  al.  /  Journal  of  Power  Sources  255  (2014)  368—376 


Table  1 

Mean  values  for  electrical  model’s  elements  for  NEW  and  the  same  cell  after  degradation  (degraded  cell). 


R0  (O) 

RS(Q) 

Cs(F) 

Ri  (O) 

Ci(F) 

New  cell 

8.2733e  -  002 

1.5115e  -  002 

8.371 6e  +  002 

3.6291  e  -  002 

4.7234e  +  003 

Degraded  cell 

1.0515e  -  001 

2.0230e  -  002 

8.371 6e  +  002 

4.2902e  -  002 

4.7234e  +  003 

Equation  (20)  is  minimized  at  each  reference  points  by  the 
optimization  method: 


(Yt  -  Voc  -  l/trans  -  Rok)  (20) 

Cell’s  terminal  voltage  and  current  are  necessary  parameters  for 
the  optimization  problem  in  equation  (20).  In  this  equation,  Vt  and 
Voc  are  associated  with  the  aged  cell  discharged  with  1C.  There  are 
three  unknown  parameters  with  two  known  SOCs.  This  equation  is 
minimized  at  both  15%  and  92%  SOC  values.  In  other  words,  there 
are  three  overall  unknown  parameters  in  two  equations.  The 
objective  function  in  two  reference  points  has  plenty  of  local 
minima  and  thus  optimization  algorithms  can  easily  get  trapped.  To 
overcome  this  problem  and  for  finding  new  model  for  aged  cell, 
brute-force  approach  is  used.  This  algorithm  is  an  enumerating 
method  in  which  all  of  the  feasible  candidate  solutions  for  the 
objective  function  are  calculated  and  the  minimum/maximum  of 
the  whole  set  is  found.  For  our  problem,  limitations  presented  in 
equation  (21)  is  used  in  optimization  method. 


•9*0,  initial  <  *o  <  1  -8*o,  initial 

•9*5,  initial  <  <  2 Rs  initial  (21) 

•9*/5  initial  <  <  2 R/  initial 

where  R0, initial,  Rs, initial  and  Ri, initial  are  the  corresponding  values  of 
R0 ,  Rs  and  R\  for  a  brand  new  cell  or  the  latest  values  of  EKF  model.  In 
the  optimization  algorithm,  £>5,  C3  and  g3  parameters  are  incre¬ 
mented  within  the  assigned  range  and  the  objective  value  is 
calculated  for  each  set  of  values.  The  set  in  which  the  mean  value  of 
error  has  the  minimum  value,  is  selected  for  placing  in  the  new 
model.  The  Cusabie  parameter  is  updated  while  the  cell  is  under  load 
using  the  method  suggested  in  Ref.  [35  .  The  estimated  SOCs  ob¬ 
tained  from  EKF  method  are  used  instead  of  the  real  SOCs. 

Based  on  the  foregoing  explanations,  this  method  can  be  sum¬ 
marized  in  the  following  steps: 


1)  Discharge  the  cell  with  constant  current  measuring  cell’s  ter¬ 
minal  voltage. 

2)  Calculate  derivative  of  the  voltage  measured  in  step  one. 

3)  Assign  SOC  92%  and  15%  as  two  reference  points  to  voltage 
derivative. 

4)  Run  optimization  algorithm  to  obtain  updated  model. 

5)  Insert  updated  model  to  IEKF. 

Fig.  5  presents  implementation  flowchart  of  the  model  adaptive 
IEKF  algorithm.  It  is  worthy  to  mention  that  similar  cells  are  filtered 
in  the  first  run.  The  range  of  resistor  values  shown  in  equation  (21 ) 
are  limited  to  the  following  bounds: 

•9^o,  initial  <  *0  <  initial 

•9*s,  initial  <  <  ll*s,  initial  (22) 

•9*/)  initial  <*/<!•!*/,  initial 


4.2.  IEKF  approach 

In  the  EKF,  the  Kalman  gain,  Kk,  represents  the  relative  impor¬ 
tance  of  the  error,  with  respect  to  prior  estimated  xk  equation 
(23).  Basically,  the  gain  controls  how  much  trust  is  there  for  the 
measurements  over  the  estimation. 

=  Xk  +  Kk(yk-Hxk)  (23) 

The  EKF  doesn’t  have  accurate  estimation  of  states  at  first  steps. 
The  Kalman  gain  guarantees  that  measurements  are  reliable  source 
of  updating  estimations.  After  200-300  steps,  when  there  is  almost 
a  good  estimation  of  the  states,  its  value  does  not  alter  significantly. 
Based  on  these  explanations,  the  IEKF  for  SOC  estimation  in  EV 
applications  is  presented.  The  linearized  system  expressed  math¬ 
ematically  based  on  equations  (5)  and  (6)  for  IEKF  is  presented  in 
Fig.  6.  This  figure  presents  a  system  with  the  estimation  of  the 
improved  extended  Kalman  filter.  In  Fig.  6,  Kk  stands  for  Kalman 
gain  and  blocks  containing  z-1  are  single  time  step  delay  boxes.  For 
estimating  states,  the  Kalman  gain’s  effect  for  first  200  steps  follows 
equation  (23).  In  the  next  steps,  since  there  is  a  good  estimation  of 
Li-ion  battery’s  states,  effect  of  residual  is  alleviated  by  reducing 
Kalman  gain’s  effect  on  equation  (23).  It  is  achieved  by  raising  the 
Kalman  gain  to  the  power  of  2.5  as  shown  in  equation  (24). 

%  =  X/7  +  (Kk- 2.5)  (yk  -  Hx, 7)  (24) 


4.3.  SOC  definition  and  estimation  of  Li-ion  cells  connected  in  series 
pack 


As  mentioned  above,  parameters  of  the  new  model  and  each 
cell’s  usable  capacity  can  be  easily  obtained  by  running  optimiza¬ 
tion  for  each  cell  in  the  pack  while  discharging  the  whole  pack  with 
a  constant  current.  Using  this  approach  and  filtering  cells  with 
more  than  ±10%  differences,  SOC  for  the  pack  is  defined  as: 


SOCpack 


Qisabl  e,i'SOQ 
Sj=i  Qisable,i 


(25) 


Fig.  4.  Voltage  derivative  for  a  new  cell  and  an  aged  cell  (discharged  with  1C). 
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Updated 
Electrical  model 
for  IEKF  approach 


Fig.  5.  Implementation  flowchart  of  the  model  adaptive  IEKF  algorithm. 


where  ns  is  the  number  of  cells  in  series  and  Cusabie,i  and  SOQ 
denote  usable  capacity  and  SOC  for  the  ith  cell,  respectively.  In  order 
to  estimate  the  SOC  of  the  pack  defined  in  equation  (25),  a  single 
unit  cell  model  for  the  pack  is  presented. 

The  single  unit  cell  model  for  the  pack  is  similar  to  the  model 
introduced  in  Fig.  1  with  following  equations  for  the  parameters  of 
battery  model: 

R0  pack  =  bi  SOC4  +  b2 SOC3  +  b3 SOC2  +  b4SOC  +  5,1  (26) 


assumed  that  those  capacitors  are  constant  and  equal  for  all  cells, 
Q.pack  and  Q  pack  are  calculated  from  equations  (29)  and  (30). 

IEKF  approach  is  used  for  SOC  estimation  of  the  pack.  Since  the 
pack’s  terminal  voltage  is  measured  at  each  time  step,  equations 
(23)  and  (24)  are  replaced  with  equations  (31)  and  (32) 
respectively. 

Equivalently,  xk  is  calculated  using  the  following  equation  ob¬ 
tained  from  IEKF: 

k  <  200  :  xk  =  xk  +  (yk  -  nsHxk )  (31) 


Vns  p; 

d  —  r  p—Cjsoc  2-w'=l  3,1 

%,pack  —  Cl"  + 

I  is 


Ri,  pack  =gie-&soc+g4SOC  + 


n< 


(27) 


(28) 


,  ,  (Kk.  2.5) 

k  >  200  :xk  =  xk+  — L  (yk  -  nsHxk )  (32) 

fLS 

Fig.  7  presents  a  relation  of  model  estimation  for  each  cell  and 
SOC  estimation  for  the  pack. 


C, 


s,pack 


Q 

tls 


(29) 


Q.pack  =  £  (30) 

where  ns  is  the  number  of  cells  in  series  and  b5  i,  c3  i  and  g3  i  are  the 
identified  values  for  the  ith  cell  of  the  pack  obtained  in  model 
adaptive  approach  as  discussed  before.  About  Cs  and  Q,  since  it  is 


5.  Experimental  results 

The  proposed  method  is  tested  in  room  temperature  on  a  bat¬ 
tery  pack  consisting  of  120  LiFeP04  cells,  connected  in  series.  These 
cells  are  A123  Systems’  APR18650ml  LiFeP04  battery  with  1.1  Ah 
nominal  capacity.  The  environment  temperature  and  thus  the 
initial  temperature  of  the  cell  is  assumed  to  be  25  °C. 

A  new  electrical  model  for  each  cell  is  found  after  applying 
proposed  optimization  algorithm  on  the  aged  battery  pack.  Next, 
SOC  is  estimated  by  application  of  IEKF.  The  SOC  estimation  for  the 


Fig.  6.  System  with  the  estimation  of  the  improved  extended  Kalman  filter  (IEKF). 
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Fig.  7.  Relation  of  model  estimation  for  each  cell  and  SOC  estimation  for  the  pack. 

battery  pack  is  compared  with  observed  SOC  which  is  based  on 
coulomb  counting  method. 

In  the  first  test,  to  show  advantage  of  model  adaptive-IEKF,  the 
battery  pack  is  charged/discharged  with  current  profile  presented 
in  Fig.  8a.  Voltage  error  estimation  and  SOC  error  estimation  of 
this  approach  and  model  adaptive-EKF  presented  in  Fig.  8b  and  c, 
respectively.  Model  adaptive-EKF  benefits  from  the  same 
approach  and  equations,  except  that  during  the  whole  test, 
equation  (31)  is  only  used  to  estimate  states  instead  of  applying 
equations  (31)  and  (32).  It  is  clear  that  model  adaptive-IEKF 
method  is  more  reliable. 

For  next  test  the  normalized  New  European  Driving  Cycle 
(NEDC)  case  study  is  used  to  evaluate  proposed  method.  This 
driving  cycle  consists  of  four  urban  driving  profiles  followed  by 
another  Extra  Urban  Driving  Cycle  (EUDC)  profile.  Based  on  equa¬ 
tion  (23),  the  pack’s  initial  SOC  is  96%.  The  current  profile  is 
consecutive  without  resting  time  till  the  pack’s  estimated  SOC 
reaches  10%.  In  Fig.  9a  the  current  profile  of  one  simulation  cycle  is 
depicted.  Fig.  9b  shows  difference  between  measured  voltage 
variation  for  the  pack  and  estimated  voltage  by  model  adaptive- 
IEKF.  SOC  estimation  error  is  presented  in  Fig.  9c.  In  this  figure, 
SOC  estimated  by  model  adaptive-IEKF  is  compared  with  observed 
SOC  calculated  by  coulomb  counting  method.  This  test  stops  as 
soon  as  estimated  SOC  for  the  pack  reaches  10%.  At  the  end  of 
experiment,  the  weakest  cell  of  the  pack  has  9.8%  SOC,  which  it  is 
still  safe  and  is  not  considered  as  over  discharged.  Mean  error  of 
SOC  estimation  is  0.4%. 

For  next  test,  a  similar  analysis  with  normalized  United  States 
Federal  Test  Procedure  (FTP-75)  is  designed.  FTP-75  is  a  more 
realistic  driving  cycle  compared  to  the  NEDC.  The  FTP-75  is  known 
as  a  transient  driving  cycle  and  represents  a  driving  in  an  urban 
environment  with  frequent  stops  [36  .  Current  profile  is  shown  in 
Fig.  10a.  Initial  SOC  for  this  test  based  on  equation  (23)  is  95%.  Test  is 


continued  with  the  same  current  profile  till  the  pack’s  estimated 
SOC  reaches  10%.  Fig.  10b  shows  voltage  estimation  error  using 
proposed  approach.  SOC  estimation  error  obtained  from  model 
adaptive-IEKF  is  depicted  in  Fig.  10c.  When  the  test  is  stopped  by 
reaching  10%  of  estimated  SOC  for  the  battery  pack,  the  weakest  cell 
in  the  branch  still  has  9%  SOC.  Mean  error  of  SOC  estimation  is  1.1%. 
Based  on  tests  results,  it  is  clear  that  presented  approach  has  more 
reliability  besides  its  accuracy.  By  reliability,  we  mean  that  none  of 
the  cells  in  the  pack  will  be  over  discharged  using  the  presented 
method  in  SOC  estimation. 

6.  Discussion  and  conclusion 

In  order  to  have  a  reliable  energy  source  and  distributing  power 
in  EVs/FiEVs,  a  battery  management  system  needs  to  have  an  ac¬ 
curate  inline  estimation  of  SOC  for  Li-ion  cells  in  the  battery  pack. 
EKF  is  an  accurate  method  for  estimating  SOC  for  Li-ion  batteries. 
This  estimation  is  accurate  while  the  model  predicts  the  cells’ 
behavior  properly.  As  battery  pack  ages,  the  model  cannot  predict 
the  battery  pack’s  output  accurately  and  EKF’s  estimation  will  not 
be  reliable  any  more.  To  overcome  this  drawback  and  improve  the 
accuracy  and  reliability  of  EKF  estimation,  a  novel  method  for 
updating  the  battery  pack’s  model  is  proposed.  The  method  used 
for  packs’  SOC  estimation  is  called  model  adaptive-IEKF.  This 
method  is  an  inline  and  on-board  method  for  updating  electrical 
model  for  a  Li-ion  battery  pack.  Based  on  experimental  results,  the 
concluding  remarks  can  be  made  below: 

(1 )  This  approach  uses  an  equivalent  electrical  battery  model  for 
describing  the  electrical  behavior  of  a  single  Li-ion  cell.  The 
model  describes  the  dynamic  characteristics  of  the  cell 
consisting  one  resistor  and  two  RC  branches  and  a  voltage 
source  related  with  SOC,  all  connected  in  series.  The  same 
equivalent  electrical  model  is  used  to  make  a  single  unit  cell 
model  which  describes  the  battery  pack’s  electrical  behavior. 

(2)  Brute-force  search  method  is  used  as  an  optimization  algo¬ 
rithm  for  finding  electrical  model  for  each  new  cell  in  the 
battery  pack.  The  same  method  is  employed  for  finding  new 
model  of  pack’s  cells  in  the  aging  process.  Mean  values  of 
identified  parameters  of  all  cells  in  the  pack  are  used  to  make 
a  single  unit  cell  model  for  the  pack.  This  model  decreases 
computational  burden  significantly. 

(3)  An  IEKF  based  algorithm  for  SOC  estimation  using  single  unit 
cell  model  is  presented  to  estimate  SOC  for  the  Li-ion  battery 


Fig.  8.  a.  Current  profile  for  first  test.  b.  Voltage  estimation  error  based  on  EKF  and  IEKF  discharged  with  current  profile  presented  in  Fig.  8a.  c.  SOC  estimation  error  based  on  EKF 
and  IEKF  discharged  with  current  profile  presented  in  Fig.  8a. 
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Fig.  9.  a.  Normalized  NEDC  current  profile,  b.  Voltage  estimation  error  based  on  IEKF  discharged  with  Fig.  9a  current  profile,  c.  SOC  estimation  error  based  on  IEKF  discharged  with 
Fig.  9a  current  profile. 


Fig.  10.  a.  Normalized  FTP-75  current  profile,  b.  Voltage  estimation  error  based  on  IEKF  discharged  Fig.  10a  current  profile,  c.  SOC  estimation  error  based  on  IEKF  discharged  Fig.  10a 
current  profile. 


pack  accurately.  This  estimator  can  predict  SOC  and  terminal 
voltage  for  Li-ion  battery  pack. 

(4)  A  Li-ion  battery  pack  with  120  cells  connected  in  series  is 
tested  under  normalized  current  profile  of  FTP-75  and  NEDC. 
Experimental  results  provide  the  pack’s  SOC  and  terminal 
voltage  error  less  than  1.5%,  which  presents  robustness  and 
reliability  of  the  proposed  approach. 

(5)  In  presented  method,  as  each  cell’s  impedance  and  usable 
capacity  in  the  battery  pack  is  estimated,  SOH  for  each  cell 
can  be  predicted  as  well.  In  sum,  this  approach  can  predict 
SOC  and  SOH  simultaneously. 
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